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A variational formulation incorporating the full Navier-Stokes equations is used to 
identify initial perturbations with finite kinetic energy Eq which generate the largest 
gain in perturbation kinetic energy (across all possible time intervals) for plane Couette 
flow. Two differen t representative flow g eom etries are chosen co r respo nding to those 
used previously bv lButler fc Farrc ll (1992) and Monokrousos et all ( 2011 ). In the former 
(smaller geometry) case as Eq increases from 0, we find an optimal which is a smooth 
nonlinear continuation of the well-known linear result at Eq = 0. At Eq = E c , however, 
completely unrelated states are uncovered which trigger turbulence and our algorithm 
consequently fails to converge. As Eq — > E+, we find good evidence that the turbulence- 
triggering initial conditions approach a 'minimal seed' which corresponds to the state 
of lowest energy on the laminar-turbulent basin boundary or 'edge'. This situation is 
repeated in the Monokrousos et al. ( 201ll) (larger) geometry albeit with one notable 
new feature - the appeara nce of a nonlinear optimal ( as found recently in pip e flow by 
Pringle fc Kerswell ( 2010l ) and boundary layer flow bv lCherubini et al. ( 2010h ) at finite 
Eq < E c which has a very different structure to the linear optimal. Again the minimal 
seed at Eq = E c does not resemble the linear or now the nonlinear optim al. Our results 
support the first of two conjectures recently posed bv lPringle et al. ( 2011 ) but contradict 
the second. Importantly, their prediction that the form of the functional optimised is not 
important for identifying E c providing heightened values are produced by turbulent flows 
is confirmed: we find the what looks to be the same E c and minima l seed using energy 
gain as opposed to total dissipation in the iMonokrousos et all (|2011l ) geometry. 



1. Introduction 

The investigation of hydrodynamic stability is one of the canonical problems of fluid 
dynamics. A particularly interesting archetypal flow is plane Couette flow (PCF), where 
the flow is between two parallel plates moving at a relative velocity 2U, separated 
by a distance 2 ft., with thus a characteristic Reynolds number Re = Uh/v where v 
is the kinematic viscosity. PCF is linearly stable for even large Re ( Romanovl (1973 )) 
yet tu rbulence has been observed experimentally as low as Re = 325 (jBottin fc Chate 
( 19981 )). It has been hypothesized that transient perturbation growth, due to the non- 
normality of the underlying linear operator of the Navier-Stok e s equations, may explai n 
this disconnect. Several authors (for cxample lGustavssonl (|199ll ): lButrer fc Farrelll (119921 ); 
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IReddv fe Henning son1 (jl993t n demonstrated that substantial transient kinetic energy gain 
G(T) = E(T)/E(0) (where E(T) is the infinitesimally small kinetic energy at the final 
time T) could be achieved by a linear (infinitesimal) optimal perturbation (LOP) (see 
ISchmi d (|2007h for a review). Proponents of such an essent i ally li near mechanism for en- 
ergy growth point to the Rcynolds-Orr equation (jSchmidl (|2007i )) as an indication that 
energy growth is a linear effect, as this equation shows that dE/dt is independent of the 
nonlinear advective terms in the Navier-Stokes equations. The argument is that if a real 
flow is seeded with a 'small', yet finite amplitude perturbation with the same structure 
as a LOP, the transient energy gain of the LOP could be sufficiently large to 'push' the 
perturbation into the 'nonlinear regime' and hence trigger transition. 

However, this line of thinking is based on a couple of implicit assumptions: that the 
'entrance' into the nonlinear regime of this perturbation will lead to turbulence; and 
that the linear optimal perturbation (LOP) is still the 'best' choice for the growth of 
nonlinear perturbations. The latter assumption can be explicitly probed by posing the 
optimal growth problem for initial perturbations of finite amplitude which can affect the 
base f low as they grow. This has been done recently fo r pipe flow ([Pringle fc Kerswell 
(|2010[) ) and boundary layer flow (|Cherubini et all (|201fj) ) with both studies discovering 
the existence of a nonlinear optimal perturbation (NLOP) which has a very different 
structure to the LOP and outgrows it beyond a small but finite energy thr eshold. 

With more of an eye on reaching turbulence, iMonokrousos et all (|201 lh posed a dif- 
ferent problem for PCF by maximising the total energy dissipation over a long but fixed 
time period. They purposely looked for a turbulent end state at the end of their opti- 
misation window and t hen worked downwa r ds in initial energy to identify the threshold 
for transition. Earlier, iPringle fe Kerswelll f)2010h had failed to identify this threshold 
by workin g upwards in i nitial energy because of convergence issues. However, a follow- 
up study (jPringle et all (|201l[ ). henceforth referred to as PWK11) with a more efficient 
code run at higher resolution succeeded in identifying a (converged) nonlinear optimal 
for larger initial energies. At Eq = Ef a n, they again found a failure to converge but 
noticed that this corresponded to their optimisation algorithm encountering turbulent 
(end) flows. The conclusion was that this failure energy Ef a u is sufficiently large to en- 
able an initial perturbation to undergo the transition to turbulence: i.e. Ef a u E c , the 
energy threshold of transition (PWK11 actually conjectured that Ef a u = E c : see Con- 
jecture 1 below). In dynamical systems parlance, this initial perturbation is then in the 
basin of attraction of the turbulence or more generally (if the turbulence is actually not 
an attractor but a chaotic saddle), has crossed the 'edge', a hypersurface w hich separates 
initial conditions wh i ch be c ome turbulent f r om th o se with relam i narise (jltano fe Toh 
( 200ll ): ISkufca et al] (|2006l ): ISchneider efaH (|2007t ): iDuguet et all (|2008l )L The picture 
then put forward by PWK11 is that as Eq increases from 0, the E = Eq hypersurface 
in phase space intersects the edge for the first time at Eq = E c and that the initial 
perturbation which corresponds to their (gcnerically unique) intersection at Eq = E c is 
the 'minimal seed' for triggering turbulence. This seed is 'minimal' in the sense that it 
is the lowest energy state on the edge and therefore represents the most energy efficient 
way of triggering turbulence by adding an infinitesimal perturbation to it. PWK11 also 
find evidence to suggest that the NLOP tends to the minimal seed associated with this 
loss of convergence (and transition to turbulence) as Eq — > E~ . They summarise their 
thinking as two conjectures. 

" Conjecture 1: For T sufficiently large, the initial energy value Ef a u at which the 
energy growth problem first fails (as Eq is increased) to have a smooth optimal solution 
will correspond exactly to E c . 

Conjecture 2: For T sufficiently large, the optimal initial condition for maximal en- 
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ergy growth at Eq = E c — e 2 converges to the minimal seed at E c as e —> 0. " 



In this paper, we wish to investigate the validity of these conjectures in the context of 
PCF. Tw o sets of geometry and R eynolds numbers are considered for PCF, one discussed 
in each of Butler fc Farrell ( 19921) (henc eforth referred to as BF92) - a relatively narrow 
spanwise domain with Re = 1000 - and iMonokrousos et al. ( 201 lh (henceforth referred 
to as Mil) who used a domain with double the spanwise extent at Re = 1500. By 
considering these different situations, we are able to investigate whether there is anything 
generic that can be said about the progression of optimal perturbations starting with the 
(infinitesimal) LOPs at Eq = 0, through NLOPs as Eq increases to the minimal seed 
at Eq = E c . Of principal interest will be whether this approach can identify E c and 
the form of the minimal seed either directly (by smooth evolution of the optimal as 
Eq — > E~) or indirectly (by failing to converge). By considering the Mil geometry but 
choosing to maximise the energy gain rather than total energy dissipation, we also assess 
the sensitivity of the procedure to the exact choice of optimising functional. 

From the technical perspective, we also take this opportunity to further develop the 
variational formulation to include optimisation over T, the duration of the observation 
window or 'target time'. This means we are then able to identify the initial perturba- 
tion which achieves the highest gain possible over all T wit h the corresponding o ptima l 
final time now an i n teresting output. All p r evious studies dPringle fe Kerswell ( 2010t ); 
Pringle et al. ( 2011 ); Cherubini et al. ( 2010() ; Monokrousos et al. ( 201ll )) chose to work 
with a pre-defined T over which to perform their optimization. Conceptually, letting T be 
an output of the optimisation seems a significant advance yet operationally, it is requires 
only a small adjustment in the algorithm. 

The plan of the paper is as follows. In section [5] we briefly present the variational 
framework and discuss how the new target time optimisation is carried out. The following 
two sections, [3] and |H discuss the results obtained for the BF92 and Mil situations 
respectively. A final section [5] then discusses the results in light of the above-quoted 
conjectures of PWK11 and draws a number of conclusions. 



2. Lagrangian framework 

We seek the initial disturbance of kinetic energy Eq to the laminar flow which attains 
the largest energy growth G(T) := E(T)/Eq a time T later while evolving under the 
Navier-Stokes equations, remaining incompressible and respecting the applied boundary 
conditions. Here E(T) := |(u(T), u(T)), with the angle brackets denoting, 

(v,u) :=IjT v tu dV, (2.1) 

where f denotes the Hermitian conjugate, and V is the volume of the domain V. The 
flow configuration considered is PCF with coordinate system such that the streamwise 
direction is x, the wall normal direction is y and the spanwise direction z. The x and z 
directions are assumed to be periodic and the separation between the walls (2h) is used 
to scale length so that their positions are given by y = ±1. The speed difference between 
the walls (2U) scales the velocity so that the non-dimensionalised background Couette 
flow is U(y) = ye x and the Reynolds number Re := Uh/u. 

The functional to be cxtrcmiscd is the energy gain which, when constrained by the 
Navier Stokes equations, the initial energy value £"(0) = Eq and incomprcssibility, pro- 
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duces the Lagrangian 

[d t u + N(u) + Vp, v] - [V.u, q] - (hu Q ,u )-E Q ) c+ (u - u(0), v ) 



TV is the nonlinear operator 



(2.2a) 



N(ui) := UjdjUi + uidiUj + UjdjUi — -prdjdjUi , (2.26) 

Re 

and square brackets denote a time average of the inner product, 



1 



T 



[v,u] :=- J (v,u)dt. (2.3) 

In the Lagrangian, v, q, Vo and c are Lagrange multipliers, Uo is the initial value of the 
perturbation velocity u and U the background Couette flow. While not strictly necessary 
to divide our cost functional by Eq we found it easier to tune our algorithm by doing so. 

Taking first variations of the Lagrangian with respect to v, q, v and c and setting 
them to zero recovers (respectively) the constraints of the Navier Stokes equations, in- 
comprcssibility, the initial kinetic energy of Eq and initial state Uo = u(0), 

— = <9 t u + iV(u) + Vp := , (2.4) 
ov 

SC , . 

— = V.u:-0, 2.5 
dq 

— = u - u(0) := 0. (2.6) 
Sj£ 1 

— = -(u ,u ) -E a :=0. (2.7) 
oc 2 

First variations with respect to the physical variables yields a complementary set of 
adjoint equations, 



^ = d t w + ATt( v , u) + Va + ( i - v) \ t=T + (v - v ) | f=0 := , (2.8) 
du \E J 

SC , . 

^ = V.v:=0, (2.9) 



5p 
SjC 
Su 



vo - cu := 0. (2.10) 



Here, 



N\vi, u) := 9j (%«j) - WjSfUj + 9j (J7j«i) - Uj9iC/j + —djdjVi (2.11) 



1 

can be identified as the adjoint of N and v, vo and q are the adjoint variables of u, Uo and 
p. Equation (|2.8[) is in reality three equations. The first part, 9 t v + N'(v, u) + Vo = 0, 
must be satisfied at all times and is the adjoint Navier Stokes equation. Since the full 
Navier-Stokes equations have been imposed, the adjoint operator depends on the velocity 
field u. The sign of the diffusion term is also reversed and therefore the adjoint equation 
can only be solved backwards in time. The second part of (|2.8j) . (u/Eq — v) \ t= T = 0, 
is a terminal condition, linking our physical and adjoint variables and needs only to be 
satisfied at time T. The third part, (v — vo) \t=o — 0, is a condition linking Vo and v(0), 
which must be satisfied at t = 0. 

Further to previous recent formulations ( Pringle fc Kerswell ( 2010t ): Pringle et al\ ( 201ll ): 
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T. The first variation with respect to T yields the simple relation 

provided u is incompressible and satisfies the Navier-Stokes equations at t = T. Our 
algorithm then proceeds as follows. We first start with a suitable guess for the optimal 
initial condition, Uo and a target time T. We then time march our initial condition to 
time T using the Navier-Stokes equations and use (u/Eq — v) \t=r = to 'initialise' 
the adjoint equations which are then solved backwards in time to calculate Vo. This 
procedure ensures that all the variational equations are satisfied apart from (|2.10p and 
(|2.12p . If the current value for Uo is optimal then (|2.10[) will be satisfied: on the other 
hand if (|2.10p is not satisfied, it provides an estimate for the gradient SjC/Svlq. Using this 
gradient we then use a method of steepest ascent to update our guess for Uo, while c 
is simultaneously calculated by ensuring that our new initial condition has an energy of 
Eq. Once a new value of Uo is obtained, T is updated by integrating the Navier-Stokes 
equations forward in time using the updated Uo as an initial condition until a maximum 
value of E{t) is reached. The time of this maximum is taken as the new value of T and 
(12. 12[) is then satisfied. 



3. BF92 geometry 

The underlying objective of this paper is to investigate how optimal initial conditions 
for energy growth some T later change as a function of Eq. One specific issue is whether 
there is always an energy range below E c where a NLOP is the optimal (a NLOP being 
an initial condition qualitatively different in structure and gain from the LOP). A second 
is whether Ef a u = E c and a third is examining the form of the optimal as Eq — > E c from 
above or below. 

Results are first presented from a geometry studied previously in the linear regime 
(Eq — > 0) in BF92: a periodic box with dimensions L x = 27t/0.49 = 13.66, L y = 2 
and L z = 2tt/1.9 = 3.31 (or 4.08tt x 2 x 1 05tt) with Re = 1000. A modified version of 
the Diablo CFD solver, (|Tavlor fc Sarkaj (|2008l )) which is spectral in x and z and finite 
difference in y, was used to solve the forward and adjoint equations using a resolution of 
128 x 256 x 32 in x, y and z respectively. For sufficiently low energies we found that the 
optimal perturbation was extremely similar to the LOP in both gain and structure and as 
result was named a 'quasi linear optimal perturbation' (QLOP). The QLOP achieved a 
maximum gain of approximately 1100 at T = 125 (in units of h/U). With increasing but 
still small Eq , the gain and optimal time of the QLOP remains fairly constant as shown 
in figure [Ha). However, beyond a certain energy threshold, (approximately 2.3 x 10~ 6 ) 
there is a sudden and large jump in the gain achievable, as shown in figure HJb) (note 
the change of ordinate scale). In addition to the much higher gain at Eq = 2.3 x 10~ 6 , 
it is noticeable that the optimal time tends to very large values as Eq approaches this 
transition energy from above. The initial conditions found by our algorithm there clearly 
evolve into a turbulent state given the much higher target-time kinetic energy values 
and the highly disordered endstate. This implies that E c < 2.3 x 10~ 6 is the threshold 
energy for transition. To examine convergence, (S£/Suq, SC/Suq) 2 /G is plotted in figure 
[2]against iteration for Eq = 5.0 x 10~ 7 < E c and Eq = 5.0 x 10~ 6 > E c . 

For the smaller initial energy case (figure (2ja)) we see that after 10 iterations the gain 
has plateaued and that the value of the normalised gradient has dropped by 10 orders of 
magnitude, which suggests that the QLOP is converging well. 
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Figure 1. (a) Gain (blue circles) and associated optimal time (red crosses) against Eq. 
E c ~ 2.3 x 10 -6 is marked by vertical green dashed line. LOP gain is the horizontal blue solid 
line, LOP optimal time is the horizontal red dashed line. (6) Gain (blue circles) and associated 
optimal time (red crosses) against Eo for Eo > 2.3 x 1CF 6 . 




Iterations Iterations 



FIGURE 2. G (blue solid line) and {SC/Suq, SC/Suq)^ /G, (red dashed line) plotted against 
iteration for (a) Eo = 5.0 x 10~ 7 (6) Eo = 5.0 x 10~ 6 . The QLOP presented in (a) appears to 
be converging well, whereas there is no convergence in (b). 



Conversely, for the higher energy case (figure EJb)) while the gain appears to plateau, 
the gradient is failing to decrease in size indicating the algorithm is not converging. 
In reality, because of the turbulent nature of the flow at time T we would not expect 
convergence to be possible, as a very small change in the initial condition is likely to 
produce a significant change in the final state. Despite this lack of convergence, the 
algorithm is successful in finding initial states which trigger turbulence when Eq > E c . 

In the picture of PWK11, one unique initial condition - the minimal seed - should 
emerge as the limiting state for turbulence-triggering initial conditions as Eq —> E^ . To 
examine the dynamical route states close to the minimal seed take to turbulent disorder, 
we have considered in detail a turbulent seed found at Eq = 2.3 x 10 -6 as it evolves 
in time. Figure EJa) plots the kinetic energy and dissipation rate against time and in 
[3fb) the same quantities are plotted for the 'rescaled' turbulent seed whose initial energy 
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Figure 3. Gain (blue solid line) and dissipation (red dashed line) against time for (a) the 
turbulent seed at Eo = 2.3 x 10 -6 and (b) rescaled turbulent seed at Eo = 2.2 x 10 -6 . Both 
initially behave similarly achieving a gain of approximately 1000. They maintain this energy 
level for an extended period of time, until t ~ 150. After this the turbulent seed has a noticeable 
spike in dissipation, indicating a transition to turbulence, whereas the rescaled turbulence seed 
decays away. 



is Eq = 2.2 x 10 -6 . The behaviours of the turbulent seed and rescaled turbulent seed 
are very similar up to t ~ 75 h/U. Beyond this time, however, the kinetic energies of 
the two initial perturbations begin to differ significantly, with the rescaled turbulent 
seed eventually decaying so the flow relaminarises whereas the turbulent seed triggers 
transition at T ~ 270 h/U. (Reassuringly, if the rescaled turbulent seed is used to initiate 
the optimizing procedure at Eq = 2.2 x 10~ 6 , the algorithm converges to the expected 
QLOP result.) 

We observe in figure [3] that both flows spend an extended period of time at an interme- 
diate (perturbation) energy level before going their separate ways. This is because both 
initial states are close to the edge (but on 'opposite sides') and spend some time tracking 
it while being gently repelled (in opposite directions). To confirm this, the turbulent 
seed and its rescaling can be used to refine the initial condition so that it stays neare r 
to the edge fo r longe r (|ltano fc Tohl (|200ll ); Iskiifca et all (|20CJ6h ; ISchneider et all (|2007l ); 
Dueuet et all (|2008l )). In fi gure [31 this refinement is carried out to track the edge for 
t = 400 h/U showing that the edge state (attr acting state for edge-co nfined dynamics) 
has constant energy (consistent with the work of lSchneider et all ( 20081 ) who treat a PCF 
system 4-7T x 2 x 2-7T albeit at Re — 400 and find a steady edge state). 

Figure \E\ shows a side by side comparison of the contours of the perturbation stream- 
wise velocity at x — at four different times for the QLOP at E = 2.2 x 10~ 6 , the 
minimal seed at E = E c (at least to the accuracy of figured]) and the turbulent seed 
(which clearly is close to the minimal seed) at E c < Eo — 2.3 x 10 -6 . These plots demon- 
strate that while the minimal seed evolves towards the edge state, a small increase in 
its initial energy will lead to transition. As Eo — > E c from above (below), the time to 
transition (relaminarisation) tends to oo due to the extra time needed to evolve upwards 
(downwards) in energy away from the edge. It is also clear that the QLOP is completely 
different from the minimal seed. 

Finally it is worth examining the 3D structure of the time-evolving QLOP, the minimal 
seed and the turbulent seed just above the edge. In figure [6] we plot iso-contours of the 
streamwise velocity for times 0, 150, 250 and 350 (in unit of h/U). From the plots it is 
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Figure 4. Perturbation energy against time for various initial states close to the edge. Upper 
bound of edge state blue, lower bound red. Every 100 time units, the edge state is rescaled to 
produce new upper and lower bounds. The minimal seed at E c (blue line) stays on the edge 
(over this time period) and gets attracted to the edge state which emerges as having constant 
energy. The turbulent seed at Eq = 2.3 x 10~ 6 is shown as the first blue/orange dashed line and 
the rescaled turbulent seed is the first red/green dashed line. 



clear that the minimal seed is initially quite localized but quickly 'unpacks' itself into 
a series of a streamwise streaks. This unpacking process appears to be achieved by the 
well-known Orr mechanism followed by the lift-up mechanism. The minimal seed flow 
then remains in this configuration, whereas the streaks destabilise and there is transition 
to turbulence for the higher energy turbulent seed 'above' the edge. 

To summarise, in this geometry and at this Re, it appears that a well-converged NLOP 
does not exist prior to our algorithm uncovering turbulence-triggering initial conditions 
at Eq = Ef a u. Our algorithm only fails to converge if there are turbulent seeds present 
so Ef a u ^ E c and we find no evidence for inequality consistent with PWKll's first 
conjecture. The minimal seed (as is apparent in figures [5] and [S]) is qualitatively different 
from the QLOP, and so the optimals do not converge to the minimal seed as Eq — > E~, 
a clear counterexample to PWKll's second conjecture. 



4. Mil geometry 

In this section we examine a second, larger geometry of dimensions 4-7T x 2 x 2tt (es- 
sentially twice as wide as that in BF92) at a higher Reynolds number Re = 1500. We 
demonstrate in this geometry that now a NLOP exists at energies below E c and inves- 
tigate whether in the limit Eq — > E~ it converges to the minimal seed. Choosing the 
geometry and Reynolds number used by Mil has the added benefit that we can com- 
pare our results to those obtained using an entirely different functional. Mil optimized 
the total dissipation over a long time interval rather than the energy gain achieved at a 
specific target time. We find a critical energy value E c = 3.3 x 10 -7 , plotted in figures [7] 
(a) and (b), which agrees well with Mil, who find 3 x 10~ 7 < E c < 4 x 10 -7 (see their 
figure 1). Note eo in Mil is Eq here as || \\e in their equation (1) is strictly a kinetic 
norm with a i included (Monokrousos, personal communication). Our calculated time 
for transition at Eq = 4.0 x 10~ 7 is approximately 200 not too dissimilar from the time 
of 150 in Mil. This suggests that the particular choice of optimizing functional is not 
important for the calculation of a minimal seed (or more accurately to lose convergence) , 
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Figure 5. Contours of streamwise velocity u at times 0, 150, 250, 350 for 
QLOP (left), minimal seed (centre) with Eq = E c and turbulent seed (right) 
for E = 2.3 x 10" 6 > E c . Contour levels are: going down the left column 
(min,spacing,max) = (-6, 2,6) x 1(T 5 , (-0.1,0.05,0.1), (-0.1,0.02,0.1) and (-0.06,0.02,0.06); 
going down the centre (min, spacing, max) = (— 2, 0.5, 1.5) x 10~ 3 and (—0.1,0.05,0.1) subse- 
quently; going down the right column (min,spacing,max) =(-2, 0.5, 1.5) x 10" 3 , (-0.1, 0.05, 0.1), 
(-0.6,0.1,0.5) and (-0.6,0.2,0.6). 



provided the functional attains heightened values for turbulent flows (as discussed in 
PWK11). 

Figure [TJ a) indicates that three energy regimes exist in this geometry rather than the 
two in BF92. As before, below a certain initial energy value, a QLOP is selected and 
above a critical energy E c initial conditions significantly different from the QLOP trigger 
turbulence. Between these two energy regions, however, there now exists a range of initial 
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Figure 6. Iso surfaces of streamwise velocity u, at 60% of maximum and minimum values, for 
the QLOP at Eo = 2.2 x 1CP 6 (left), the minimal seed (centre) and a turbulent seed above the 
edge at E = 2.3 x 1CT 6 (right), at times 0, 150 ,250, 350. 



energies where our algorithm generates an initial condition different from the QLOP (in 
the sense that it appears to have a qualitatively different spatial structure) - see figured! 
Using the nomenclature described in the introduction, we call th is qualitatively different 
optimal perturbation a NLOP after IPringle k Kerswelll (|2010h and PWK11. Figure M 



contrasts the convergence for the NLOP with the non-convergence in the turbulent seed 
region E n > E c . Only five points are plotted in figure |H] as for additional iterations 
(however small we made our step size in the direction of the gradient) it was not possible 
to find a new u with a gain that improved on the previous iteration. 

As a consequence of the kinetic energy gains of the NLOP and QLOP becoming very 
similar around £o = 1.0 x 10~ 7 , the cross-over between the NLOP and QLOP is hard to 
pinpoint. If the NLOP is used to initialise the algorithm for energies slightly above E c , 
the algorithm is found to converge to an initial condition very similar to the NLOP. In 
fact, the algorithm started with random noise will still sometimes converge to the NLOP 
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Figure 7. (a) Gain, E{T)/E(Q), of QLOP (blue circles) and NLOP (orange circles) and as- 
sociated optimal time of QLOP (red crosses) and NLOP (green crosses), T, against Eq. E c is 
marked by vertical green dashed line. LOP gain horizontal blue solid line, LOP optimal time 
horizontal red dashed line. (6) Gain, E(T)/E(0), (blue circles) and associated optimal time (red 
crosses), T, against Eq. E c is marked by vertical green dashed line. 




Figure 8. Iso surfaces of streamwise velocity u, at 60% of maximum and minimum values, for 
(a) QLOP at E = 5.0 x 10~ 8 and (6) NLOP at E = 3.2 x 10~ 7 . It is clear that the NLOP is 
distinct from the QLOP. 

for values of Eq twice as large as E c . As a consequence, approaching E c from above 
proved a better strategy. Random noise was used at Eq ~ 2.5 x 10~ 6 rj 7.5E c to find a 
turbulent seed and then this was used sequentially to initiate the algorithm as Eq was 
gradually decreased. This experience clearly emphasizes the main hazard of nonlinear 
optimisation: it is easy to get stuck near local maxima. Although not a cure, an obvious 
strategy to reduce this possibility is to look for robustness of result over a suite of initial 
conditions. 

Figure [7] also indicates that the turbulent seeds show the same trend, as in the BF92 
geometry, with regards to the optimal target time, namely that it increases drastically as 
Eq —> E^ . Wc conclude that the turbulent seeds remain near the edge for an even greater 
period of time than the turbulent seeds in the BF92 geometry. This may be because they 
are closer to the edge and/or that the edge is less repelling. As before, we have traced 
the edge up to t = 400 h/U using a slightly rescaled turbulent seed to find a similar 
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Figure 9. G (blue solid line) and (SC/Suo, 5£/5uq) 2 /G (red dashed line) plotted against iter- 
ation for (left) E = 1.0 x 10~ 7 and (right) E = 3.3 x 10~ 7 . The NLOP (left) appears to be 
converging well, whereas the calculation which throws up turbulent seeds does not. 



plot to figure 2] (not shown) . A comparison of cross sectional streamwise velocities for 
the NLOP at E — 3.2 x 10 -7 , the minimal seed and the turbulent seed at 3.3 x 10~ 7 
in figure [TTJ again emphasizes their different temporal evolutions despite being so close 
energetically. Also interestingly, the NLOP is localised in the cross-stream direction and 
is not dissimilar from the minimal seed although they are clearly not the same. This is 
made obvious by comparing their streamwise structure: see the top row of figure [TO 

Figure [TT] also shows how the NLOP unpacks into a series of streamwise streaks. 
An examination of the early time suggests that it is a combination of the Orr and 
lift-up mechanisms (as discussed in PWK11) that is responsible for the localized flow 
unpacking into streamwise streaks. The minimal and turbulent seeds also unpack in 
the streamwise and cross-stream direction producing streamwise streaks which are still 
spanwise localised. If there is sufficient energy in these streaks they are unstable (the 
turbulent seed) otherwise not (the minimal seed). By comparing the isosurfacc of the 
QLOP at time zero (figure [8]) and the isosurfaccs depicting the time evolution of the 
NLOP it is clear that the NLOP evolves into a structure which is very similar to the 
QLOP. This suggests that while at early time there exists a distinct NLOP structure, 
which is able to extract enhanced gain from the base flow by 'unpacking', it later exploits 
the same 'lift-up' mechanism as the QLOP at intermediate times. 

It is significant that the minimal and turbulent seeds are spanwisc-localised (at least 
until the turbulence is reached) in this 27r wide geometry and not in the 1.057T wide 
geometry of BF92. Of course the higher Re must be a cont r ibutor y factor but the 
spanwise dimension does seem important. IPringle &: Kerswell ( 2010t ) originally found 
an azimuthally-localised (and radially-localised) NLOP in a short pipe where one could 
talk about a 'spanwise' (azimuthal) lengthscale of 27T (radii or half-channel heights). The 
emergence of a NLOP in the wider geometry is also noteworthy. The LOP, and by def- 
inition QLOP, are global periodic states which are largely insensitive to the geometry, 
whereas the gathering evidence is that the NLOP is an attempt by the fluid to localise 
in order to maximise the energy gain for given global kinetic energy. As a result, the 
general trend should be for the energy cross-over from QLOP to NLOP to decrease with 
increasing domain size. Clearly this cross-over is above E c for the BF92 geometry at 
Re = 1000 and below E c for the Mil geometry at Re = 1500. 
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Figure 10. Contours of streamwise velocity u for NLOP Eo = 3.2 x 1CT 7 (left), approxi- 
mated minimal seed and turbulent seed at Eg — 3.3 x 10~ 7 (right), at times 0, 75, 150, 250 
and 400. At intermediate times, the minimal seed flow on the edge and that initiated from 
the turbulent seed remain relatively similar. Contour levels are: going down the left column 
(min,spacing,max)=(-10,2,8) x 10~ 4 , (-0.08,0.02,0.04), (-0.1,0.02,0.08), (-0.1,0.02,0.08), 
and (—0.06,0.02,0.06); going down the centre (min, spacing, max) = ( — 10, 2, 8) x 10~ 4 , 
(-0.1,0.05,0.1), (-0.15,0.05,0.1), (-0.2,0.05,0.05) and (-0.3,0.05,0.1); going down the 
right column (min,spacing,max) = (-10, 2, 8) x 10" 4 , (-0.1,0.05,0.1), (-0.15,0.05,0.1), 
(-0.3,0.05,0.1) and (-0.8,0.2,0.8). 



The evolution shown in the right column of figure [TT] looks very similar to that shown 
in figure 4 of Mil (albeit at a slightly different initial energy) indicating that we have 
generated an approximation to the same (presumably unique) minimal seed. This is 
further supported by the aforementioned correspondence in their and our estimates for 
E c . Beyond validating each others results (which is important for nonlinear optimisation 
problems), this points to an insensitivity in the choice of the functional to be maximised 
for finding the minimal seed. There is one proviso, of course, that the functional must 
be selected so that it detects turbulent flows by assuming large values as discussed in 
PWK11. 
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Figure 11. Iso surfaces of streamwise velocity u, at 60% of maximum and minimum value, 
for NLOP Eo — 3.2 x 10 -7 (left), approximated minimal seed (centre) and turbulent seed at 
Eq — 3.3 x 10~ 7 (right), at times 0, 75, 150, 250 and 400. The minimal and turbulent seeds are 
initially localized but quickly unpack into streamwise streaks, which are stable for the minimal 
seed but unstable for the turbulent seed ultimately leading to breakdown. 
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In this paper, we have sought the disturbance to plane Couette flow of a given finite ki- 
netic energy E which will experience the largest subsequent energy gain G = E(T)/E 
where the time of maximum gain T is an output of our variational formulation. Two 
geometry- Re situations have been considered: (L x x L y x L z , Re) = (4.087T x 2 x 



1.057T , 1000) as it was used for the original calculations of linear optimals in lButler fc Farrell 



( 1992h and (4ir x 2 x 2ir, 1500) for which analogous calculations hav e recently been per 



forme d optimising the total dissipation over a specified period in iMonokrousos et 
(|201lh . Our results can be summarised as follows. 



(a) A nonlinear optimal (NLOP) distinct from the 'nonlinearised' linear optimal (QLOP) 
exists only in the wider geometry at the Re considered. 

(b) In both situations, there exists an energy Ef a u beyond which the variational algo- 
rithm no longer converges due to the existence of turbulence-triggering initial conditions. 
This means Ef a u ^ E c , the energy above which turbulence can be triggered. PWKll's 
first conjecture is that Ef a u = E c if the energy hypersurface is sufficiently sampled and 
we find nothing to contradict this. 

(c) The QLOP or NLOP are not found to converge to the minimal seed (the distur- 
bance of lowest energy which can trigger turbulence) as Eq —> E~ contradicting PWKll's 
second conjecture. 



(d) The failure of our variational algorithm to optimise energy gain appears to give 
the same estimat e for E c and the minimal s eed as optimising the total dissipation over a 
long time period IMonokrousos et al\ (|20 1 lh . This confirms the 'robustness of failure' of 
the variational approach discussed by PWK11 providing the functional to be optimised 
assumes large values for turbulent flows. 



The underlying motivation for all these types of optimising calculations is the hope 
that the 'optimal' perturbation at Eq < E c for an appropriately selected functional 
bears some relation to disturbances of lowest energy which actually trigger turbulence. 
By optimising the energy gain, PWK11 found this was at least approximately the case for 
their pipe flow set-up. Here though, this is clearly not true in either PCF situation studied. 
Ironically, however, the variational procedure does identify these turbulence-triggering 
disturbances and the critical energy for transition E c but only indirectly by failing to 
converge. Importantly, we have also found evidence that this revealing failure to converge 
is insensitive to the exact functional selected providing the functional takes on heightened 
values for turbulent flows as argued in PWK11. Taken together, the variational approach 
adopted here seems to offer a fairly robust new theoretical tool to examine the nonlinear 
stability of fluid flows. 

Many applications suggest themselves, but here we just note one - assessing the stabi- 
lizing or destabilizing influence of applied flow perturbations or controls. Normally, this 
would be attempted either by investigating the linearised operator around the base (lami- 
nar) flow or by carrying our exhaustive numerical simulations. The current work suggests 
a third way where the movement (in phase space) of the laminar-turbulent boundary to- 
wards (destabilisation) or away (stabilisation) from the base flow is investigated. We hope 
soon to report on some calculations along these lines. 
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